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ABSTRACT 

I demonstrate by means of high resolution cosmological simulations, which include 
modelling of a two-phase interstellar medium, that the dominant mechanism for trans- 
porting heavy elements from the proto-galaxies into the IGM is the merger mechanism 
as discovered by Gnedin & Ostriker. Direct ejection of the interstellar gas by super- 
novae plays only a minor role in transporting metals into the IGM: for a realistic 
cosmological scenario only a small fraction of all metals in the IGM is delivered by 
the supernova-driven winds, while most of all metals in the IGM are transported by 
the merger mechanism. As the result, the metallicity distribution in the IGM is highly 
inhomogeneous, in agreement with studies of the QSO metal absorption systems, and 
the predicted metallicity distribution of Lyman-alpha absorbers as a function of their 
column density is in excellent agreement with the observational data. 

Key words: cosmology: theory - cosmology: large-scale structure of universe - galax- 
ies: formation - galaxies: intergalactic medium ~ intergalactic medium: evolution - 
hydrodynamics - methods: numerical - quasars: absorption lines 



1 INTRODUCTION 



It is now well established that the Lyman-alpha forest ab- 
sorbers, thought to mostly consist of density peaks in the 
low density intergalactic medium away from (proto-) galaxies 
(Cen et al. 1994; Zhang et al. 1995; Petitjean, Mucket, & 
Kates 1995; Hernquist et al. 1996; Miralda-Escude et al. 
1996; Wadsley & Bond 1996; Bi & Davidsen 1996; Hui, 
Gnedin, & Zhang 1996), contain a measurable, about 1% 
of solar, amount of heavy elements ( "metals" ) at a redshift 
as high as three and even beyond. An early indication to this 
extent by Meyer & York (1987) is now confirmed quantita- 
tively by new high-resolution Keck data (Womble, Sargent, 
& Lyons 1995; Songaila & Cowie 1996), who found that as 
a column density of the Lyman-alpha absorbers decreases, 
smaller and smaller fraction of all absorbers contain a mea- 
surable amount of heavy elements, indicating highly inho- 
mogeneous distribution of metallicities in the intergalactic 
medium (IGM). What is the origin of those metals and their 
inhomogeneous distribution? 

Within the boundaries of standard astrophysics, there 
is only one "factory" that produces metals - stars. Stars 
lose mass via winds and explode as supernovae, polluting 
the interstellar medium with metals. However, from a cos- 
mological point of view, there is a long way from the in- 
terstellar medium, which has a cosmological overdensity 5 
(5 = p/p — 1) of about lO"* -7- 10^, to the intergalactic 
medium with overdensities of about 10 and below; any mech- 
anism that transports metals from the proto-galaxies into 



the IGM should therefore be powerful enough to overcome 
the gravitational force associated with enormous overdensi- 
ties of (proto-)galaxies.Q 

For a long time since the pioneering works by Couchman 
& Rees (1986) and Dekel & Silk (1986; see also Miralda- 
Escude & Rees 1997 for a recent review) it has been thought 
that the supernovae blowing out gas from the Population III 
objects would do the job. However, recent numerical three- 
dimensional simulations (Ostriker & Gnedin 1996; Gnedin 
& Ostriker 1996) as well as semi-analytical results (Haiman, 
Rees, & Loeb 1996) convincingly demonstrated that Pop HI 
objects, formed via molecular hydrogen cooling, are so few 
(mostly because the molecular hydrogen is so fragile), that 
they are unable to produce the metallicity in excess of about 
10"^ solar. 

Nevertheless, supernovae can still disrupt normal. Pop 
II, proto-galaxies and pollute the IGM with metals, as ar- 
gued, for example, by Miralda-Escude & Rees (1997). If 
it were the only mechanism for transporting metals from 
the deep potential wells of proto-galaxies into the diffuse 
IGM, it would be the final answer to the question: "What 
is the metal transport mechanism in the IGM?". However, 
as was discovered by Gnedin & Ostriker (1996), there exists 



I avoid using the term "galaxies" , which has a well established 
meaning, since I will discuss properties of self-gravitating objects 
at a range of redshifts, including high redshifts z ~ 10, when no 
galaxies existed. 



© 0000 RAS 



2 Gnedin 



another metal transport mechanism, which is based solely 
on gravitational hierarchical clustering. Namely, when two 
proto-galaxies merge with a sufficiently large relative veloc- 
ity in a close to head-on collision, dark matter halos and stel- 
lar components of two halos penetrate through each other, 
whereas gas components of both objects collide in a gigan- 
tic shock. The pressure force axising from this interaction 
is comparable to the gravitation force from the daxk matter 
and baryons of two objects. Therefore, the collisionless com- 
ponent (the dark matter and stars) and the coUisional com- 
ponent (the interstellar gas) experience forces that signifi- 
cantly differ in their directions and amplitudes, and, as the 
result, both components may have different spatial distribu- 
tions in the process of merging. Because the gas constitutes 
only a small fraction of the total mass of a proto-galaxy, it 
often happens during those violent mergers that some frac- 
tion of the gas (typically a few percent), displaced relative 
to the dark matter, becomes gravitationally unbound and 
gets ejected from the proto-galaxy into the IGM. This gas, 
which was a part of the interstellar medium in the past, 
and is heavily enriched with metals, gets mixed up with the 
surrounding IGM, polluting it with heavy elements. As the 
result, the metal distribution in the IGM becomes highly 
inhomogeneous: while regions around the proto-galax;ies are 
enriched with heavy elements, the centers of large voids, well 
removed from the places where merging occurs, stiU contain 
the pristine material (an example of how this mechanism 
works is shown in Gnedin & Ostriker 1996 and will not be 
presented in this paper). 

One can imagine other ways to transport metals from 
the proto-galaxies to the IGM. For example, a fraction of 
young stars can be lost during a merger, and going off later 
as supernovae, may enrich the IGM. I consider this process a 
part of merging mechanism, and will not distinguish it from 
the case when the metal enriched gas rather than stars was 
lost in the merger. After all, since stars are subject to the 
same (purely gravitational) force as the dark matter is, and 
the gas is not, it is much easier to unbound a substantial 
fraction of the total gas component of two merging proto- 
galaxies rather than of their stellar component. 

The whole process of the metal enrichment of the IGM 
begins in stars: stars make metals in the first place in the 
gas with the cosmic overdensity of J ~ 10^"; blowing up 
as supernovae, they enrich the ISM {d ~ 10®) with metals; 
metals from the ISM get transported into the IGM [5 ^ 
10) cither by the very same supernovae or by the merging 
mechanism as explained above. It is therefore the purpose 
of this paper to investigate the relative importance of those 
two metal transport mechanisms within the framework of a 
reasonable cosmological model. 

Since the transition between the ISM and the IGM is 
somewhat uncertain, and in order to avoid introducing an- 
other, "intermediate" , component for the cosmic gas, I will 
call "the ISM" the gas that is gravitationally bound to a 
particular proto-galaxy, and "the IGM" the rest of the cos- 
mic gas (i.e. gas, which is not bound to a particular self- 
gravitating object). 

Several attempts have been made before to include ef- 
fects of supernovae into cosmological hydrodynamic simu- 
lations (c.f. Cen & Ostriker 1992; Navarro & White 1993; 
Katz, Weinberg, & Hernquist 1996). However, all previ- 
ous works attempted to model the ISM as a sinpte-phase 



medium, simply injecting thermal and kinetic energy of su- 
pernovae into the cosmic gas. This approach is strongly de- 
pendent on the numerical resolution of a simulation. For a 
low resolution simulation, injection of thermal energy may 
produce a significant effect. However, as the resolution of 
a simulation increases, and higher and higher density re- 
gions are being resolved, the cooling time in the gas becomes 
smaller and smaller, and any amount of the thermal energy 
that is directly injected into the cosmic gas, is radiated away 
before it has had a change to influence the dynamics of the 
gas. In reality however, the thermal energy gets injected in 
supernova remnant bubbles (the second phase in the gas), 
that are too hot and do not cool as efficiently. 

On the other hand, the kinetic energy directly injected 
into the gas can propagate a large distance from the source 
of injection, because the unsufficient resolution of a simula- 
tion prevents one to resolve the snow-plow phase of the su- 
pernova remnant evolution (Cioffi, McKee, & Bertschinger 
1988; McKee 1990), during which most of the kinetic en- 
ergy of the supernova ejecta is absorbed. Thus, the effect of 
the finite numerical resolution can lead to both underesti- 
mating and overestimating the effect of the supernovae on 
the dynamics of the ISM, depending on the resolution of a 
simulation and details of the implementation of the super- 
nova feedback. Because the resolution of a high-resolution 
Lagrangian simulation is not uniform across the simulation 
volume, and varies from point to point, it is extremely diffi- 
cult to disentangle the physical cfi^cct of supernovae feedback 
from numerical artifacts produced by the finite resolution of 
a simulation. 

It is therefore logical to attempt to model the multi- 
phase ISM in order to account properly for the supernovae 
feedback. Again, the finite resolution of existing cosmolog- 
ical hydrodynamic simulations does not allow to resolve 
properly the evolution of a single supernova remnant, but 
the multi-phase ISM can be incorporated in a cosmological 
simulation as "sub-cell" physics, similar to the way the star 
formation is included in simulations. In other words, for a 
sufficiently large region, which contains many supernovae, it 
is possible to account for the total effect, a population of su- 
pernova remnants have (on average) on the dynamics of the 
interstellar gas (energy injection, extra pressure, reduction 
in the star formation rate etc). 

Simulations reported on in this paper attempt to do 
precisely that. The main difference between this work and 
previous investigations is that a too-phase medium is now 
allowed for in a simulation, i.e. at each point in a simula- 
tion the cosmic gas has two temperatures (the "cold" phase 
outside the supernova remnants and "hot" phase inside the 
supernova remnants) and two densities (which can be ex- 
pressed as the total density and the volume fraction occu- 
pied by the hot phase). The current implementation explic- 
itly assumes however that the population of supernova rem- 
nants is moving together with the rest of the cosmic gas, i.e. 
the two phases have the same velocity on a large (~ 1 kpc) 
scale. The velocity itself is however affected by the feedback 
from the supernova remnants. The specific recipe for incor- 
porating physics of the two-phase ISM into a cosmological 
hydrodynamic code is presented in Appendix A. 

Finally, a few words about interpreting results of a cos- 
mological simulation. Any "sub-cell" , phenomenological de- 
scription of a physical process like supernova feedback or 
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star formation inevitably introduces new parameters that 
describe the physical process under consideration. It is cur- 
rently unrealistic (from the point of view of computational 
resources) to perform a full-scale extensive parameter study, 
and thus it should be kept in mind that any cosmological 
hydrodynamic simulation which includes "sub-cell" physics 
(including simulations reported on in this paper) should be 
considered as a plausible possibility rather than the final an- 
swer. There always remains a chance that the nature picked 
up different values of parameters than those used in a sim- 
ulation, however plausible the latter seem to be. 



2 METHOD 

I have adopted a CDM-I-A cosmological model as a frame- 
work for these investigations. The qualitative results of my 
simulation are applicable to the whole family of the CDM- 
type models, and I emphasize those qualitative conclusions 
together with specific quantitative results for the model I 
use. The cosmological parameters are fixed as follows: 



h 



0.35, 
0.65, 
0.70, 
0.055, 



which is one of the cosmological models adopted by the 
Grand Challenge Cosmological Consortium (GC^) as the 
"best-bet" models {http://arcturus.mit.edu/GC3/IC). This 



model is similar to the model investigated by Gnedin & 
Ostriker (1996) except that Qb has a higher value, in con- 
cordance with the low primordial abundance as inferred by 
Tytler & Buries (1996). 

Since the merging mechanism can only be studied by a 
numerical simulation, I use the SLH-P'^M cosmological hy- 
drodynamic code as described by Gnedin (1995), Gnedin 
(1996), Gnedin & Bertschinger (1996), and Gnedin & Os- 
triker (1996). However, since the goal of this paper is to 
study the effects of supernovae on the metal transport in the 
IGM, a proper numerical treatment of the physical effect of 
supernovae has to be incorporated in the code. Appendix 
A explains how the physics of the multi-phase interstellar 
medium can be incorporated into a cosmological hydrody- 
namic code on scales large compared to the size of one super- 
nova bubble. Three parameters control the evolution of the 
ISM in this regime: -B51, the average energy of a supernova 
measured in units of 10^^ erg, /imf, a factor parametriz- 
ing the excess (lack) of metal production compared to the 
Miller- Scalo IMF, and t;tcrm, the termination velocity of the 
supernova shock at which the interior of the bubble is ion- 
ized. In order to minimally span the parameter space and un- 
derstand how supernovae effect the IGM on different scales, 
six simulations have been performed. Each simulation has 
64^ dark matter particles, 64^ gas cells, and a number of 
stellar particles which are being formed during a simulation. 
The softening parameter for all simulations is set to 1/10, 
corresponding to the dynamical range of 640. In addition to 
the new ingredient, the physics of the supernova, physical 
modelling in all runs includes all components described in 
Gnedin (1996) and Gnedin & Ostriker (1996), to which pa- 
pers the reader is referred for more details about physical 



Table 1 


. Numerical Parameters 








Run 


Box size 


Spatial res. 


£^51 


/imp 


crm 


Al 


l/i-l Mpc 


l.e/i^^kpc 





4 


NA 


Bl 


Ih-i Mpc 


1.6/i~ikpc 


1 


4 


cs 


CI 


Mpc 


l.eft-^kpc 


1 


4 


30km/s 


A3 


3h-^ Mpc 


4.7/1-1 kpc 





4 


NA 


B3 


3h'^ Mpc 


4.7/1^1 kpc 


1 


4 


cs 


C3 


3h'^ Mpc 


4.7/1-1 kpc 


1 


4 


30km/s 



modelling. Here I only note that I allow for the fully three- 
dimensional treatment for the evolution of the dark matter 
and the two-phase cosmic gas, as well as for the formation 
and destruction of molecular hydrogen and all other stan- 
dard physical processes for a gas of primeval composition, 
following in detail the ionization and recombination of all 
species in the ambient radiation field. The spatially and fre- 
quency dependent radiation field, in turn, allows for sources 
of radiation (quasars and massive stars), sinks (due to con- 
tinuum opacities) , self-shielding of the dense gas lumps, and 
cosmological effects. In regions which are cooling and col- 
lapsing (beyond the numerical resolution of the code) I al- 
low the formation of point-like "stellar" subunits, permit- 
ting them to release radiation and explode as supernovae. 
Supernovae inject thermal and kinetic energy in the ISM 
(via the second gaseous phase - the supernova remnants), 
and pollute the surrounding gas with metals, which are in 
turn considered in the treatment of cooling. 

Parameters of all runs are summarized in Table ^ In 
order to span a large enough range in mass scales, two sets of 
runs are performed: the first set with the lh~^ Mpc box size 
and the total mass resolution of lO^'^fe-'^M©, and the second 
set with the 3h~^ Mpc box size and the total mass resolu- 
tion of IO^ '^/i-^Mq . The resolution in baryons is a factor 
of 0.055/0.35 = 0.16 finer, and for runs with the Ih'^ Mpc 
box size, it is W^'^ h~^MQ, which is about the Jeans mass at 
z ~ 15 for this model. Those runs, therefore, account for all 
initial fluctuations that are present in the baryons, i.e. they 
do not miss any small scale power in the initial conditions 
(but those simulations do miss the small scale power that is 
generated by the nonlinear dynamics at later times on the 
scales below the resolution limit). 

Simulations within each set have had exactly the same 
initial conditions; so that evolution of individual objects can 
be followed under different physical assumptions. Because 
it makes no sense to continue a simulation of a finite size 
beyond the redshift when perturbations with the wavelength 
of the order of the box size become nonlinear, runs with 
3/1" 1 Mpc box size are stopped at the redshift z = 4, and 
runs with Ih^^ Mpc box size are not continued beyond z — 
9. 

In each of two sets three different runs have been per- 
formed: 

run A, when no physics of the multiphase ISM (as de- 
scribed in Appendix A) has been taken into account (man- 
ifested by the parameter E51 set to zero); this run does not 
include effects of supernovae, and therefore only the merg- 
ing mechanism for metal transport operates in this case; this 
run is then used for comparison with two other runs where 
effects of supernovae are included, to assess the relative im- 
portance of two metal transport mechanisms. 

run B, which includes physics of the multiphase ISM; in 
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this run it is assumed that the supernova shock is terminated 
at the sound speed of the ambient medium. 

run C, which includes physics of the multiphase ISM; in 
this run the supernova shock is terminated at the velocity 
of 30 km/ s. 

The latter choice is somewhat arbitrary; it approximately 
maximizes the effect of supernova on heating of the low 
metallicity ISM. As has been shown by ShuU & McKee 
(1979), the supernova shock in the ISM with solar metallicity 
fails to ionize the gas at shock velocities below 90km/s. It 
is expected that for low metallicity gas at high redshift this 
value is somewhat lower (ShuU & Silk 1981), and I somewhat 
arbitrary fix it to be 30 km/ s, near the value which maxi- 
mizes the effect of supernovae. This choice also allows for the 
possibility that the star formation occurs inside the dense 
supernova swept-up shells, which are considered belonging 
to the cold phase for the shock velocities below 30km/s. 

Finally, I choose a factor /imf = 4, which implies that I 
assume a top-heavy IMF, with a factor of four higher metal 
production then in the solar neighborhood. As argued by 
Silk (1996), this value is supported by metal abundances in 
the clusters of galaxies if one assumes that all metals came 
from cluster galaxies alone. This number is an extreme, in a 
sense that if the clusters of galaxies contain additional stars 
which are distributed over the cluster (and which would be 
virtually unobservable), then /imp should be less than 4. 
Both runs CI and C3 should therefore be considered as lim- 
iting cases when all possible parameters are pushed to their 
extreme values in order to maximize the effect of supernovae. 



3 RESULTS 

3.1 The Metallicity of the IGM 

As the cosmic structure develops under the influence of grav- 
ity, baryons condense into bound objects, cool and even- 
tually form stars, which end their lives as supernovae and 
pollute the environment (i.e. ISM and IGM) with metals. 
Figure |l] shows in the upper panel the comoving star for- 
mation rate in units of Mq/ yr/{h^^ Mpc)^ for all six runs: 
the bold lines indicate runs A3-C3, and the thin lines cor- 
respond to runs Al-Cl. The solid lines mark runs A, the 
dotted lines mark runs B, and the dashed lines mark runs 
C. In all figures below line markings are the same unless 
specified otherwise. Note that the star formation declines at 
late epochs due to the final size of a simulation box (and, as 
the result, lack of large-scale power to drive the gravitational 
collapse further). 

The lower panel of the same figure shows the mass- 
weighted average metallicity of the cosmic gas (i.e. ISM and 
IGM together) as a function of redshift. I would like to note 
here, that simulations with difl^erent box sizes are reasonably 
consistent with each other in a sense that runs with the 
3h~^ Mpc box size (A3-G3) have similar star formation rates 
and metallicities at z = 9, when runs with the l/i~^Mpc 
box size (Al-Cl) become unreliable because of the missing 
large-scale power. 

As can be seen, both the star formation rates and the 
average metallicities for runs with the lh~^ Mpc box size 
are somewhat suppressed in cases when supernova effects 
were taken into account (runs Bl and CI, thin dotted and 




20 10 8 6 4 

1 + z 



Figure 1. Upper panel: the comoving star formation rate in units 
of Mq/ yr:/{h~^ Mpc)^ as a function of redshift for six runs: Al 
{thin solid line), Bl [thin dotted line), CI [thin dashed line), A3 
{hold solid line), B3 {bold dotted line), and C3 {bold dashed line). 
Lower panel: the mass weighted mean metallicity of the cosmic 
gas as a function of redshift; line markings are as on the upper 
panel. 



dashed lines) as compared to the case when no supernova 
effects were included (run Al, the thin solid line). It is the 
expected result, since as the volume fraction occupied by 
supernova bubbles increases, a smaller and smaller fraction 
of the gas is available for star formation, because the gas 
inside the bubbles is hot and is not able to cool fast enough. 
However, even on the smallest resolved scales this effect is 
not dramatic: in the most favored case (CI) both the star 
formation rate and the average metallicity are decreased by 
2: = 9 by a factor of 2.5, but not by orders of magnitude. 

Before I proceed further with more detailed compari- 
son between the two main mechanisms of metal transport 
in the IGM, it is instructive to discuss the spatial distribu- 
tion of the metallicity. Figure |^ shows the mass- (the upper 
panel) and the volume-weighted (the lower panel) proba- 
bility distributions for a fluid element to have given values 
of the metallicity and the gas density (as measured by the 
overdensity S = p/p — 1) from run Alat2: = 4(a similar 
figure is presented in Gnedin & Ostriker 1996). The imme- 
diate conclusion is that the distribution is quite broad, with 
regions at about the cosmic mean density having values of 
metallicity all the way from zero to about solar (at z = 4 \). 
The distribution narrows down in the high density limit, be- 
cause the high density gas resides inside the proto-galaxies 
and is enriched with metals. This high-density, metal-rich 
gas occupies only a small fraction of the volume, so, if there 
is a substantial fraction of the gas with metallicities as high 
as tenth of solar (at z = 4 !), most of the volume is occu- 
pied by the low metallicity gas. This wide distribution of 
metallicities is in agreement (at least qualitatively) with the 
recent conclusions by Ranch, Haehnelt, & Steinmetz (1996) 
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Figure 2. The mass- [upper panel) and volume- {lower panel) 
weighted probability distribution to find a fluid element with 
given values of the metallicity and the cosmic gas density for 
run A3 at z = 4. The bold line shows the contour corresponding 
to the probability P{Z = {Z),5 = 0)/4, solid contours mark the 
probability higher and the dashed contours show the probability 
lower than the one marked with the bold contour; contour spacing 
is logarithmic with the increment of 1/3 dex. 



and by Hellsten et al. (1997), who found that the observed 
scatter in the abundance ratios of the QSO metal absorption 
systems is inconsistent with the assumption of the uniform 
metaUicity. The required scatter in the metallicity at a given 
value of the cosmic density is about an order of magnitude, 
which is consistent with Fig. |^ for S ^ 10. 

I emphasize here that Fig. |^ presents the total cosmic 
gas, including both the ISM and the IGM. The gas with 
overdensities in excess of 10^ -r- 10'* is apparently located 
inside the proto-galaxies, and is identified with the ISM. The 
low density gas, S ^ 10, is clearly a part of the IGM (and 
is, therefore, partly pristine), with the gas at intermediate 
densities belonging to the ISM or the IGM, depending on 
whether it is bound to a particular self-gravitating object or 
not. 

In order to compare the metallicity-density relationship 
for different runs, I consider how the average metallicity 
changes with the gas density. For of all fluid elements hav- 
ing the overdensity 5 within the range from S to 5 + dS, I 
compute their metallicity Z{5) as a function of 5 and plot it 
in Figure ^ in units of the mass- weighted average metallicity 
for all six runs (in the other words, the solid bold line in the 
left panel is Fig. |^ collapsed along the vertical direction, and 
other lines are produced from the metallicity-density distri- 
butions for other respective runs; line markings are again 
as in Fig. |^). The left panel shows runs (A3-C3) at 2: = 4 
(since runs Al-Cl were stopped at z = 9), and the right 
panel shows all six runs at z = 9. 

The fact that all three lines almost coincide in the left 
panel of Fig. demonstrates that at intermediate redshifts 




1+5 1+6 

Figure 3. The average metallicity of the g as as a function of the 
cosmic gas density at z = 4 [left panel) and z = 9{right panel). 
Line markings as in Fig. ^ except for z = 4 case where no lines for 
runs Al-Cl are shown since those runs were stopped at a higher 
redshift. 



{z ~ 4) and large scales (~ 100 kpc+1 Mpc) supernovae play 
no role in transporting metals into the IGM. It is the merg- 
ing mechanism that enriches the IGM with heavy elements 
and produces a highly inhomogeneous distribution of metal- 
licities. Even at higher redshifts, and on smaller scales, su- 
pernovae play a minor role as a metal transport mechanism 
in the IGM. However, in the right panel of Fig. ^ (z = 9) one 
can notice that in run CI (the most favorable for the super- 
novae effects) there is a fraction of gas sitting at very low 
densities (about 0.01 of the mean) and having the metal- 
licity 40 times above the average (about 10% solar). This 
would be a signature of the supernova ejected gas, which 
would be highly enriched with metals and would have low 
densities, since its temperature is higher (about 10^ K) that 
the average temperature of the IGM at this high redshift 
(10^ + 10* K). However, except for this very low density gas, 
differences in the metallicities between runs with supernovae 
(B,C) and runs without supernovae (A) are not significant 
for the IGM with the density higher than the cosmic mean. 

The difference between runs Al (no supernovae) and 
CI (most efficient supernovae) is further demonstrated in 
Figure ^, which shows the mass- weighted metallicity-density 
distributions for those two runs in a fashion similar to the 
upper panel of Fig. ^ except that two different line markings 
now correspond to two runs: the solid lines show contours 
for runs Al, and the dashed lines show contours for run CI. 
The two distributions are very similar except for a small 
amount of the very low density - high metallicity gas which 
is present in run CI but is absent in run Al. This gas is 
blown out of a proto-galaxy by the supernova activity; its 
amount is insignificant: it occupies 0.6% of the total volume 
and amounts to 0.03% of the total gas mass in the universe. 
This gas is not necessarily the only gas that was ejected from 
the proto-galaxies by supernovae: there are slight differences 
in the metallicity-density distributions in runs Al and CI 
at higher densities, which are due to differences in the star 
formation rates and the supernova activity (but those two 
can not be meaningfully separated). However, in the whole, 
the supernova mechanism plays a minor role in transporting 
metals into the IGM, and by z ~ 4 gas disruption in violent 
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Figure 4. The mass-weighted probability distribution to find a 
fluid element with given values of the metallicity and the cosmic 
gas density for runs Al (solid lines) and CI {dashed lines) at 
2 = 9. The contour spacing is as in Fig. ^. The low-density - 
high-metallicity region of the plane in run CI (supernova blown 
gas) accounts for 0.03% of the total gas mass and 0.6% of the 
total volume. 



mergers is responsible for delivering of more than 99% of all 
heavy elements in the IGM. 

One may ask a question: what is fraction of the enriched 
gas in the colliding high redshift objects which gets unbound 
and ejected into the IGM? It might not be appropriate to 
talk about a single value, but it is easy to come up with 
an upper limit. Let me consider an oversimplified picture of 
the universe, which consists of high-density bound objects 
with the average metallicity Zh, which contain /h of the 
total mass. The rest of the universe is in the IGM, with 
the average metallicity zl- I also assume that bound objects 
occupy a negligible fraction of the volume. In this case, the 
volume weighted averaged metallicity will be just Zl- For 
the run A3 at z = 4, this value is = 6 x 10~^ Zq. The 
mass weighted average metallicity in this case is {Z)m = 
Zh fn + ZlJl- For run A3 at z — 4 the respective values are 
/jj — 0.2 and {Z)m = 5 x 10~^ Zq, which implies that the 
average metallicity of bound objects Zh ~ 0.23Zq. Since all 
metals in the IGM come from the gas that gets unbound in 
mergers, the fraction f^^b of this gas is 

ZlJl 



/unb — 



= 10%. 



Zh/h + Zl/l 

This number is however an upper limit, since I have assumed 
that the bound objects occupy no volume at all, because any 
clustering in the IGM will decrease this bound as well, and 
because the limited resolution of my simulations underesti- 
mates the mass fraction fn locked in bound objects. 

More than ten years ago Dekel & Silk (1986) anal- 
ysed the effect of supernovae on the gas ejection from high 
redshift proto-galaxies using the analytical approach. They 
found that supernovae can efficiently expel gas from low 



mass proto-galaxies provided that the star formation is ex- 
tremely rapid and is happening on a free-fall time scale. 
I point out here that there are two important physical ef- 
fects, which are included in the simulations presented here, 
but omitted in the Dekel fc Silk (1986) calculations. First, 
the negative feedback [eq. (A8)] limits the star formation 
rate when the fraction of the volume occupied by supernova 
remnants becomes large. This effects is ignored in Dekel & 
Silk (1986) paper. Second, Dekel & Silk (1986) ignore cool- 
ing of the interstellar gas, which leads to a substantial loss 
of energy injected by the supernova. For example, if all of 
the supernova energy was simply injected into the ISM gas 
without proper account for the two-phase ISM and super- 
nova remnant evolution, all energy would be radiated away, 
producing negligible net effect. 

One can therefore ask a question whether the present 
work confiicts with the conclusions of Dekel & Silk (1986). 
I would like to emphasize here that this paper does not pre- 
tend to claim that the supernovae do not eject gas from 
bound objects, and as I show in the following section, they 
sometimes do it very efficiently. Thus, the supernova ejec- 
tion is a mechanism for metal transport into the IGM. How- 
ever, as I explained above, there is another metal transport 
mechanism, the merger mechanism, and the simulations pre- 
sented in this paper support the conclusion that the merger 
mechanism is more efficient. The major physical reason for 
that is that the merger mechanism works irrespectively of 
the binding energy of an object and of the star formation 
rate in the object, whereas supernovae are able to disrupt 
the gas from a bound object only for objects with sufficiently 
low masses and high star formation rates. More than that, 
in order to transfer metal enriched gas into the IGM, it is 
required to expel the gas not just beyond the virial radius 
of the dark matter halo, but beyond the turn-around ra- 
dius; otherwise, the expelled gas will fall back into the dark 
matter halo in a fraction of a Hubble time. 

Therefore, it has to be kept in mind that the results 
presented in this paper do not necessarily contradict the 
conclusions of Dekel & Silk (1986), but a more careful com- 
parison will have to be performed in order to verify whether 
the simulations actually confirm the conclusions of Dekel & 
Silk (1986). 

3.2 Supernovae in the Proto-Galaxies 

As I have demonstrated in the previous section, the 
supernova-driven winds are less efficient in transporting 
metals in the IGM than the merger mechanism. Neverthe- 
less, supernova disruption of proto-galaxies does happen, as 
the difference between runs Al and CI is due to the gas dis- 
rupted by the supernovae from the proto-galaxies. In order 
to understand this process, I show in Figure ^ properties of 
all well-resolved (containing more than 100 particles of all 
kinds) bound objects identified at 2: = 9 in runs Al (the 
left column), Bl (the middle column), and CI (the right 
column). Three rows show the gas fraction, the stellar frac- 
tion, and the fraction of the hot gas within the supernova 
remnants for each object versus its total mass (the lower 
left panel is empty because the supernova effects are not 
included in run Al). The fraction of the hot gas never ex- 
ceeds about 30% since it is averaged over the whole object, 
including the outer regions where there is no star formation 
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Figure 5. Properties of all bound objects containing more than 
100 particles identified at 2; = 9 in runs Al {left column), Bl 
{middle column), and CI {right column). Shown in three rows as 
a function of the total mass of an object are the gas fraction for 
each object {upper row), the stellar fraction [middle row), and 
the fraction of the hot gas within the supernova remnants {lower 
row). The only object whose gas was completely disrupted by the 
supernova activity in run CI, and the same objects in runs Al 
and Bl, are marked with a circle in the upper row and in the 
right column. 



(and thus the supernova activity); in most of the objects it 
reaches close to 100% in the center, where star formation 
takes place. In run CI fractions of the hot gas are smaller 
that the respective fractions in run Bl, because the super- 
nova remnants are assumed to terminate earlier in run CI. 
However, the total effect of supernovae is actually larger in 
run CI because each supernova remnant excerts much larger 
pressure on the ambient ISM and ejects much more thermal 
energy than a respective supernova remnant in run Bl. 

In all three columns distributions are similar, except 
that in the third column (run CI) there is a single object 
whose gas fraction is very small, and whose gas was appar- 
ently expelled by supernovae. This object and respective ob- 
jects in runs Al and Bl are marked by a circle (since all three 
runs Al-Cl started with the same initial conditions, there is 
a one-to-one correspondence between the bound objects in 
all three runs). One can also notice at least three objects in 
the same panel lying somewhat below the main distribution. 
Those objects have also lost some of their gas due to the su- 
pernova activity, but their gas fraction is reduced by about a 
factor of 2, and not by two orders of magnitude. This would 
imply that the total disruption of a bound object by super- 
novae is a relatively rare event, but the supernova-driven 
winds are much more common. 

The evolution of the disrupted object (the bold lines) is 
presented in Figure ^ together with the evolution of the most 
massive object (the thin lines) in a simulation. Four panels 
show the total mass (the upper left panel), the stellar mass 




Figure 6. Evolutionary tracks for objects from runs Al and CI 
marked with circles in Fig. ^ {bold lines) together with evolu- 
tionary tracks for the most massive objects from the same runs 
{thin lines). Line marking is as before: solid lines show run Al, 
and dashed lines show runs CI. The gas disruption occurred at 
z = 11. Four panels show the total mass {upper left panel), the 
stellar mass {lower left panel), the gas mass {upper right panel), 
and the average gas temperature {lower right panel) as a function 
of redshift. 



(the lower left panel), the gas mass (the upper right panel), 
and the average gas temperature (the lower right panel) as 
a function of redshift. The disrupted object is shown by the 
bold dashed line. The act of disruption, which is manifested 
by a sharp drop in the gas mass and an equally sharp rise 
in the temperature occurred at z = 11. It was preceded by 
a steep rise in the object's stellar mass, which is a signature 
of the burst of star formation. For comparison, the most 
massive object (the thin lines) has similar behavior of both 
the total and the stellar mass until the moment of disrup- 
tion, with two exceptions: its stellar mass is rising slightly 
less steeply just before z=ll, and its total mass continues 
to rise, whereas the total mass of the disrupted object rose 
insignificantly from 2: = 12.5toz=ll. Itis therefore plau- 
sible to assume that the disruption of the gas contents of the 
object is a result of a rare coincidence of two events happen- 
ing almost simultaneously: a burst of star formation and a 
sharp decline in the accretion rate. It is not very likely that 
those two processes happen simultaneously, and, therefore, 
the total disruption of the gas component of a proto-galaxy 
is a rare event. The more quantitative discussion of the pro- 
cess of gas disruption is however beyond the framework of 
this paper. 

Figures serve to illustrate the process of gas dis- 
ruption. These figures show a neighborhood of the object 
marked with a circle in Fig. ^. The scale of each panel is 
200/1^^ kpc X 200/i~^ kpc in comoving units. The upper pan- 
els show the gas temperature, and the lower panels show the 
gas density for the same object in runs Al (the left panels) 
and CI (the right panels). Fig. ^ shows the object at the 
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Figure 7. A slice of the gas density (lower row) and the gas 
temperature (upper row) of comoving dimensions 200/i~^ kpc X 
200/i~^ kpc X 15h~^ kpc around the disrupted object at 2 = 11 
for runs Al (left column) and CI (right column). Stars (i.e. SLH 
stellar particles) are shown by white symbols. The small square 
at the center mark the mean cosmic density for the density panels 
and T = 10'^ K for the temperature panels. The color version of 
this plot is available online as explained in the Appendix B. 



moment of disruption, 2 = 11, and Fig. |8| shows the same 
comoving volume at z = 9. In order to illustrate the further 
evolution of this objects, runs Al and CI are continued until 
z — 7, with Fig. y showing the same comoving volume at 
z = 7. Stars are shown on the same plots with white stellar 
symbols. 

The process of disruption started at a = 11 with the 
ejection of the gas in the "north-east" direction (assuming 
that north is on the top of the figure), because the high- 
density filaments prevented the shock wave to expand into 
the north and east directions. The low density - high tem- 
perature region formed inside the shock wave. At a later 
time, 2 = 9 (Fig. ^), the shock wave overcame the confin- 
ing pressure of filaments and expanded in all directions; a 
very low gas density region formed around the proto-galaxy, 
which now contains almost entirely stars and the dark mat- 
ter (not shown). The expanding shock wave snowplowed 
a large amount of external material, and secondary shock 
waves were propagating in all directions inside the expand- 
ing bubble. As it expanded and snowplowed the surrounding 
gas, the shock wave also managed to disrupt a small fraction 
of stars from the object; stars interact with the other matter 
only gravitationally, but since the gas constitutes some 15% 
of the total mass of the object (for a considered cosmologi- 
cal model), the mass of the expanding shell is not negligible, 
and some stars were disrupted from the proto-galaxy by the 
gravitational pull of the expanding shell. By 2 = 7 (Fig. |^, 
the expanding shell dissolved and mixed with the surround- 
ing IGM, heating the gas to much higher temperature than 
it would otherwise had. 




Figure 8. The same as Fig. M except at z = 9. 
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Figure 9. The same as Fig. Q except at z = 7. The upper right 
panel is of somewhat poorer quality than other panels because of 
incomplete rebinning. 



For comparison, the left columns (run Al without su- 
pernovae) show nothing even close to this spectacular activ- 
ity. 

An MPEG video of this process as well as color versions 
of Fig. are available online as explained in Appendix B. 

3.3 Metal Abundances in the Lyman-alpha Forest 

While both the theoretical modelling and the observations 
of the metallicity of the Lyman-alpha Forest are still in the 
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Figure 10. The fraction of all peaks in the gas density which 
have the metallicity above the limiting value = 10~^Zq as 

a function of the cosmic gas density for runs Al {solid lines), Bl 
{dotted lines), and CI {dashed lines) as extrapolated to z = 3. 
The upper scale shows the average column density of a Lyman- 
alpha forest line arising from a respective density peak (Hui et al. 
1996). Points with errorbars are observational data from Womblc 
et al. (1995) and Songaila & Cowie (1996) (the data point from 
Womble et al. (1995) is corrected by 17% for Zum = W'^ Zq 
from their value of Zn^ = 10~^'^Zq using the dependence of the 
fraction of all peaks above Zii^ from run Al). 



developing stages, it is still instructive to make a rough 
comparison of the results of simulations with the observa- 
tional data from Womble et al. (1995) and Songaila & Cowie 
(1996). As was demonstrated by Hui et al. (1996), most of 
the Lyman-alpha lines with column densities below about 
10^^ cm~^ are arising from the peaks in the line-of-sight dis- 
tribution of the cosmic gas density. Figure |lo| shows by the 
bold lines the fraction of all one-dimensional density peaks 
(i.e. peaks along any one of three directions) with the metal- 
licities in excess of 1% solar as a function of the peak density 
at z — 3. Three almost coinciding lines correspond to three 
runs: A3 (the solid line), B3 (the dotted line), and C3 (the 
dashed line). Since all three runs were stopped at z = 4, I 
have extrapolated the metallicity distributions from those 
three runs to z = 3. 

Since the density of a peak strongly correlates with the 
column density of a Lyman-alpha line, it is possible to estab- 
lish a one-to-one correspondence between the peak density 
and the column density of the absorption line. Respective 
column densities are presented on the upper scale in Fig. |l^ 
for a considered cosmological model (the correspondence be- 
tween the peak density and the column density of a Lyman- 
alpha line depends on cosmological parameters). Solid dots 
with errorbars are data from Womble et al. (1995) (the 
leftmost point) and from Songaila & Cowie (1996). Since 
Womble et al. (1995) quote results for the limiting metallic- 
ity of 10~^'^Zq, I have corrected their value for the limiting 
metallicity of 1% solar using the dependence of f{Z > Znm) 
on Zuni as derived from run A3 (by 17% percent). A some- 
what arbitrary 25% errorbar was added to each data point. 

One can note that the agreement between the observa- 
tions and the simulations is pretty good, taking into account 
that no fitting has been performed to achieve this agree- 
ment. However, this agreement should only be considered as 
encouraging, but by no means definite agreement between 



the theoretical predictions and the observations. Since sev- 
eral theoretical uncertainties such as limited resolution of 
my simulations, scatter in the column density - peak den- 
sity relation etc, are not included in the analysis, it remains 
to be seen if more accurate simulations (including generat- 
ing synthetic Lyman-alpha spectra and analyzing them by 
fitting Voigt profiles) will confirm this agreement. 

It also remains to be seen how sensitive this test is. 
Since the merging mechanism is the dominant metal trans- 
port mechanism in the IGM, and the amount of merging 
strongly depends on the amplitude of density fluctuations, 
one may expect that cosmological models with different am- 
plitudes of density fluctuations on small scales (several hun- 
dred kpc) will have substantially different predictions for the 
fraction of Lyman-alpha lines above a certain metallicity as 
a function of a column density. 



4 CONCLUSIONS 

1 have demonstrated by means of high resolution cosmologi- 
cal simulations that the dominant mechanism for transport- 
ing heavy elements from the high density cosmic gas (ISM) 
into the low density cosmic gas (IGM) is the merger mech- 
anism as discovered by Gnedin & Ostriker (1996). Direct 
ejection of the interstellar gas by supernovae plays only a 
minor role in transporting metals into the IGM: for a realis- 
tic cosmological scenario (i.e. the one that has a reasonable 
amount of power on all relevant scales from the COBE scale 
all the way to the scale of about 100 kpc probed by the low 
column density Lyman-alpha forest), and at the most favor- 
able conditions, only a small fraction of all metals is directly 
expelled from proto-galaxies by supernovae, while most of 
all metals in the IGM are delivered by the merger mecha- 
nism. It does not mean that supernovae do not expel gas 
from proto-galaxies, and examples in Section 3.2 show that 
they indeed do. However, the main conclusion of this paper 
is that the other, merger mechanism is more important in 
transporting metals into the IGM (in particularly, because 
it is not a subject to the negative feedback process), and ig- 
noring it in computing the metal transport rate in the IGM 
is inadequate. 

As the result, the metallicity distribution in the IGM 
is highly inhomogeneous, in agreement with conclusions by 
Ranch et al. (1996) and Hellsten et al. (1996). 

The metallicity distribution of Lyman-alpha absorbers 
as a function of their column density is in excellent agree- 
ment (subject to the unknown uncertainty due to the limi- 
tations of theoretical modelling) with the observational data 
from Womble et al. (1995) and Songaila & Cowie (1996). 

Apparently, the effect of the supernovae depends on the 
star formation rate. One may therefore that the star forma- 
tion rate can be underestimated in the simulations, thus 
leading to a erroneous conclusion that supernovae are not 
efficient in expelling metals from the ISM. While details of 
the specific implementation of star formation in the simu- 
lations reported on here can be found in Gnedin (1996), I 
note here that the algorithm adopted in the simulation as- 
sumes that the star formation occurs on the dynamical time 
in the rapidly cooling gas at a rate, which, if being contin- 
ued until 2 = 0, would substantially overpredict the stellar 
mass density in the universe at the current epoch. Thus, the 
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star formation rate in the simulations, as all other relevant 
parameters, is pushed to the reasonable extreme in order 
to maximize the effect of supernovae. And still, even in the 
most favorable conditions, supernovae do not contribute sig- 
nificantly to the metal transport rate into the IGM at high 
redshift (z <, 4). 

This conclusion is only relevant in the range of red- 
shifts simulated in this paper, and for the class of cosmo- 
logical models whose properties are similar to the LCDM 
model considered here. The current work makes no pre- 
dictions about the relative importance of supernova driven 
winds at lower redshifts, where higher star formation rate 
and lower density of the cosmic gas in the virialised objects 
may lead to enhanced role of supernova driven outflows. 

As a side note, I mention here that since the merg- 
ing mechanism works irrespectively of the binding energy 
of a collapsed object, merging in clusters of galaxies at the 
current epoch also ejects metals in the IGM. Thus, I can 
make a prediction that, if the conclusions of this paper are 
indeed correct, the mean metallicity of the IGM at the cur- 
rent epoch can easily be as high as several percent, or even 
10% of solar 
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anonymous referee. This work was supported by the UC 
Berkeley grant 1-443839-07427. Simulations were performed 
on the NCSA Power Challenge Array under the grant AST- 
960015N and on the NCSA Origin2000 mini-supercomputer 
under the grant AST-970006N. 



APPENDIX A: INCORPORATING THE 
MULTI-PHASE INTERSTELLAR MEDIUM 
INTO A COSMOLOGICAL SIMULATION 

As has been shown by CiofB et al. (1988), a supernova rem- 
nant radius Rs in the pressure driven snowplow stage (when 
radiative losses set in) is given by the following expression: 



Rs = 25.7 pc 



p31/98 
^51 



f5/98„18/49 3/7' 

? rig (JjQQ 



(Al) 



where £51 is the energy of the supernova remnant in units 
of 10^^ erg, no is the ambient hydrogen density in units of 
lcm~^, wioo is the shock velocity vs in units of 100 km/ s, 
and 



as shown by CiofB & ShuU (1991). This radius is reached at 
the time ts after the explosion. 



t Assuming that 10% of the gas lost in mergers of clusters of 
galaxies, which contain the intracluster gas with the 1/3 solar 
metallicity, leads to the conclusion that a 3% solar metallicity in 
the IGM is produced just from clusters of galaxies at the current 
epoch. 



ts = 7.57 X 10* ys 



p31/98 



S "0 ''lOO 



(A2) 



For the very low density ambient medium no < ricrit, when 
radiative there are no losses, where 



n„it = 2.8 X 10"=* cm"^ ,„ , 

^3/2^1/2 

the supernova remnant radius is given by the following ex- 
pression: 

^1/3 

^s = 27.8pc^^. 



(A3) 



The volume fraction /hot of the ISM occupied by the 
supernova remnants is determined by the porosity Q, 

/hot = 1 - exp(-Q), 

where (McKee 1990; Silk 1996) 

1 dpt 



Q = l>- 



(A4) 



msN dt 

Here V is the 4-volume of the supernova remnant in the 
space-time, dp,/dt is the rate of star formation, and msN 
is the mass in stars formed per supernova. The 4-volume of 
the supernova remnant is given by the following expression: 



4-7r „o 
= q^RstS: 



(A5) 



where a dimensionless parameter q depends on the time- 
evolution of the supernova remnant. 

I now assume that the remnant evolution terminates at 
the shock velocity utorm- If vterm ~ l3cs, where cs is the 
sound speed of the ambient ISM, and /3 is a parameter of 
order unity (Cioffi et al. 1988), then the 4-volume of the su- 
pernova remnant until the moment of maximum expansion 
is (McKee 1990) 



1^1 



47r 

-Rsts, 



where rj = 0.3 is the power-law index of the expansion law, 
Rs ^ ig- Comparing to equation ( |A5| ), I get qi = 1/{3t]+1). 
However, in this case one must also account for the contrac- 
tion phase, which gives §2 ~ l/(4»7) (McKee 1990). The total 
q in this case (which I will call Case A hereafter) is 



qA = 



+ 



Sri + l 477 

The volume of the ISM occupied by the supernova remnants 
contains hot ionized gas, and, therefore, no stars can form 
inside the supernova bubbles. However, as was shown by 
ShuU & McKee (1979), a supernova remnant shock fails to 
ionize the ambient medium with the solar metallicity at the 
shock velocities below about 90km/s. In this case the inter- 
stellar gas snowplowed by the supernova shock with a lower 
velocity will not be ionized and may form stars even if it 
located inside the supernova remnant. In order to account 
for this possibility, I consider vterm as a free parameter, not 
necessarily close to the sound speed of the ambient ISM. 
Obviously, in no case can Vterm be smaller than Pes, but it 
can be larger. In the latter case (Case B) the contraction 
phase is not counted as being a part of hot ISM, and 

1 
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Supernova remnants affect the ambient medium in two 
ways: they exert an addition pressure on the ambient gas, 
and, when they terminate, they deposit the rest of their 
mechanical energy that was not cooled off, into the ambient 
ISM. The pressure jump at the shock wave is 

Po 7+1 ^ ' ' 

where 7 is the polytropic index for the gas, and M is the 
Mach number of the shock. 



M = 



Vtc 



CS 



For the power-law expansion law the shock velocity changes 
as t''"^, and the average relative excess pressure a supernova 
remnant exerts on the ambient ISM over its lifetime is 



Po 



7 + : 



7 + 1 <? 



where dimensionless parameter p is different for cases A and 
B: 

1 1 

Case A: pA ~ 



5r? - 1 3?7 + 1 ' 



Case B: pb 



5ri 



Pes 



3v + l 



Note, that in the latter case pb (and, therefore, the excess 
pressure) can be very large if t;tcrm ^ /3cs. 

After the supernova remnant dies, the rest of its en- 
ergy is deposited in the ISM as thermal heating. The total 
thermal energy density deposited in the ISM by a supernova 
remnant is then 



AU 

Uo 



27 
7 + 1 



2 + (7- 



(A7) 



where Uo is the ene rg y de nsity in the ambient medium. 

Equations (A4), (A6), and (A7) can now be incorpo- 
rated in the SLH-P^M code. Let the subscript denotes 
quantities that would be computed by the code in the ab- 
sence of supernova, i.e. when /hot ~ 0. The hydrodynamic 
gas pressure is corrected as P = Po + AP/hot, the star for- 
mation rate is reduced by a factor of /hot, 



dp, 
IF 



(A8) 



and the thermal energy is ejected in the gas at the rate 

fd-U^ _.rrr 1 dm, 

\~J7 1 - ^^/liot TT- 

\ at / SN msN at 

I would like to point out here that the reduction in the star 
formation rate itself depends on the porosity of the IGM, 
which itself depends on the star formation rate. This mutual 
dependence works as a powerful feedback mechanism. 

The last undefined parameter in this treatment is msN, 
the total mass in stars formed per supernova. For a solar 
neighborhood (or Miller-Scalo) IMF this value is msN = 
25OM0. I introduce a new parameter, /imf, defined as 



/iMF 



250M, 







msN 



to account for possible deviations from the Miller-Scalo IMF. 
As argued by Silk (1996), a value of /imp ~ 4 is required in 
order to account for all metals in the clusters of galaxies. 



Finally, since pB (and, therefore, AP) increases with 
increasing the termination velocity 

^tcrm, one may imag- 
ine that the effect of supernovae may be increased indefi- 
nitely by increasing Utorm. However, since both Rs and ts 
decrease with increasing Utorm, the porosity of the ISM de- 



creases with increasing iiterm as Q 



-19/7 



(or 



-24/7 



for no < ?icrit). Therefore, for a sufficiently large fterm, when 



the porosity is small, /hot = Q 



-19/7 



and the total cor- 



rection to the gas pressure (and to the thermal energy) de- 
creases with increasing vtam as AP ~ 



-5/7 



There is there- 



fore a value for iitcrm when the effects of supernova on the 
ambient ISM is maximal. 



APPENDIX B: ATTACHMENTS 

Additional illustrations to this paper are available online at 
the URL: 



http://aBtron.berkeley.edu/~gnedin/metigra.htraj 
They include: 



but with the slice width of 



• color versions of Fig. [7 j 

• color versions of Fig. 
Ih^^ kpc instead of 15h^^ kpc; 

an MPEG video of the evolution of the object discussed 



in §3.2 



a color image showing the merging mechanism at work. 
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